Universal phase diagram of a strongly interacting Fermi gas with unbalanced spin 

populations. 
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We present a theoretical interpretation of a recent experiment presented in ref. [l4| on the 
density profile of Fermi gases with unbalanced spin populations. We show that in the regime of 
strong interaction, the boundaries of the three phases observed in [3] can be characterized by two 
dimensionless numbers rj a and r)p. Using a combination of a variational treatment and a study of 
the experimental results, we infer rather precise bounds for these two parameters. 

PACS numbers: 03.75.Hh, 03.75.Ss 
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INTRODUCTION 

In fermionic systems, superfluidity arises from the pair- 
ing of two particles with opposite spin states, a scenario 
first pointed out by Bardeen Cooper and Schrieffer (BCS) 
to explain the onset of superconductivity in metals. For 
this mechanism to be efficient, the Fermi surfaces asso- 
ciated with each spin component need be matched, and 
soon after the seminal BCS work the question of the ef- 
fect of a population imbalance between the two states was 
raised. At the time, it was understood that pairing and 
superfluidity could sustain a certain amount of mismatch, 
above which the system would undergo a quantum phase 
transition towards a normal state [l|. The original work 
of Fulde, Ferrel, Larkin and Ovchinnikov, who proposed 
the existence of Cooper pairing at finite momentum, was 
later generalized to trapped systems [2| . Alternative sce- 
narios were also proposed, including deformed Fermi sur- 
faces interior gap superfluidity phase separation 
between a normal and a superfluid state through a first 
order phase transition [B[ , BCS quasi-particle interactions 
or onset of p-wave pairing @. When the strength of 
the interactions is varied, a complicated phase diagram 
mixing several of these scenarios is expected [jj. 

However, due to the absence of experimental evidence, 
these scenarios could never be tested experimentally un- 
til the subject was revived by the possibility of reach- 
ing superfluidity in ultra-cold fermionic gaseous systems 



lfj,lll[. Contrarily to usual condensed matter systems, 



spin relaxation is very weak in cold atoms, and this allows 
one to keep spin polarized samples for long times. This 
unique possibility led to the first experimental studies 
of imbalanced Fermi gases at MIT and Rice University 
12\ lil [13 ] . These results triggered a host of theoretical 
work aiming at explaining the various results observed 
by the two groups 15, lil ]. 



One remarquable feature of ref. [14| is the observation 
of three different phases in the cloud. At the center, the 
authors observe a superfluid core, where the densities of 
the two spin states are equal, then an intermediate nor- 
mal shell where the two states coexist and finally an outer 



rim of the majority component. In the present paper, 
we show that, though performed in a trap, the observa- 
tions of MIT can offer valuable information on the phase 
diagram of a strongly interacting Fermi gas with unbal- 
anced populations. In a first part we will present a brief 
overview of the simplest free space scenario for transition 
from a paired superfluid to a pure normal state, through 
a mixed phase. Focusing on the disappearance of the 
minority component, we will present a variational study 
of the problem of a single minority particle embedded in 
the Fermi sea of majority atoms. Finally, we will show 
that the comparison with experiments allows for a rather 
precise determination of the transition thresholds. One 
of the key point is that, contrarily to previous works, we 
rely on universal thermodynamics [T3| as well as "exact" 
experimental or Monte Carlo results, without the need 
of the mean field BCS ansatz often used in other publi- 
cations, an approach similar to that of [Hj]. 



HOMOGENEOUS SYSTEM 

Before addressing the case of trapped fermions, let us 
first discuss their free space phase diagram. We consider 
zero temperature fermions of mass m with two internal 
states labelled 1 and 2. Within a quantization volume V 
and in the limit of short range interactions, we can write 
the hamiltonian of the system as 
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Here, e*, = h 2 k 2 /2m, afc j(T is the annihilation opera- 
tor of a species a particle with momentum fc, and g\, is 
the bare coupling constant characterizing inter-particle 
interactions. It is related to the s-wave scattering length 
of the system a by the Lippmann-Schwinger equation 
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We note that only interactions between particles of oppo- 
site spins are taken into account, due to the Pauli prin- 
ciple which forbids s-wave scattering of atoms with iden- 
tical spin. In this paper, we assume we are working at 
the unitary limit where \a\ — oo and, using the grand 
canonical ensemble, we wish to find the ground state of 
the grand-potential 5 = H — fiiNi — Here Ni 

is the particle number operator for species i and /_*, is 
the associated chemical potential (we take species 1 as 
majority hence fi\ > /j,?). The grand potential can be 
expressed as the function of the volume V and pressure 
P of the ensemble according to (H) = — PV. In other 
words, searching the ground state of the system is equiv- 
alent to searching the phase with the highest pressure 
P. 

To start our analysis, we note first that two exact 
eigenstates of S can be found quite easily. First, when 
the gas is fully polarized, we recover the case of an ideal 
Fermi gas for which we know that 

1 /2m\ 3/2 5/2 

Second, let us now consider the exact ground state 
SF) M of the balanced grand potential S' = H — fi(N± + 
N2), describing a superfluid with chemical potential \i. 
This potential commutes with the number operators, 
hence the ground state can be searched as an eigenstate 
for both Ni and N 2 , with iVi|SF) M = 7V 2 |SF} M . We 
check readily that |SF)p, with p = (/xi + M2)/2, is still 
an eigenstate of the unbalanced S, by noting that S = 
S' + (fix — H2){Ni — iV 2 )/2. At unitarity the pressure of a 

balanced Fermi gas reads Ps = 2 (2m/£/z 2 ) 3 ^ 2 /1 5 / 2 /I5n 2 , 
where £ ~ 0.42 is a universal parameter whose de- 
termination has attracted interest of both theoreticians 
[S O IS, HH and experimentalists [HI, . In the case 
of mismatched chemical potentials, the pressure of this 
fully paired superfluid state is therefore 

The evolution of Pjy and Ps is presented in Fig. [1] 
as a function of 77 = ^2/^1- We see that they cross for 
i] c = (2£) 3 / 5 — 1 ~ -0.10, marking the instability of 
the superfluid against large population imbalances [la ]. 
However, since we only compare the energy of the fully 
paired state to the one of the fully polarized ideal gas, 
the real breakdown of superfluidity could very well hap- 
pen for some r\ larger than rj c . We know this is actually 
the case, since in ref. , the authors observed an inter- 
mediate normal phase, containing atoms of both species. 
From universality at unitarity, the phase transition from 
the fully paired to the intermediate phase, and then from 
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FIG. 1: Comparison of the pressure of the various phases, 
normalized by the pressure Pn of the fully polarized ideal 
Fermi s function of the chemical potential mismatch 

7] = nz/fii. Dotted line: fully polarized phase. Full line: fully 
paired superfluid phase. The fully paired and fully polarized 
states meet for 7/ c ~ —0.099. Dashed line: sketch of the 
intermediate normal phase. r) a and 77/3 designate the universal 
chemical potential thresholds for this phase. 



the intermediate to the fully polarized normal phase are 
given by conditions /J2/M1 = Vai and ^2/^1 = Vi3> where 
rj a and rip are two universal parameters we would like to 
determine as precisely as possible. 

Noting that the transition from the fully paired state to 
the intermediate one must happen before the transition 
to the fully polarized phase, we see graphically that we 
have necessarily r\ a > r\ Cl and, similarly, r/p < rj c . The 
upper bound on rjp can be further improved by noting 
that at the threshold between the normal mixture and 
the fully polarized ideal gas, there are only a few atoms 
of the minority species. In principle, the value of r/p 
should then be found by studying the N+l body problem 
of a Fermi sea of N particles 1, in presence of a single 
minority atom. To address this problem, we use here a 
variational method inspired from first order perturbation 
theory, where we expect the ground state of the system 
to take the form 



-0 o |FS) + ^0 fc , q |fc,q), 

k.q 

where |SF) is a non interacting majority Fermi sea plus a 
minority atom with momentum, and \k,q) is the per- 
turbed Fermi sea with a majority atom with momentum 
q (with q lower than Uf) excited to momentum k (with 
k > fop). To satisfy momentum conservation, the minor- 
ity atom acquires a momentum q k. The energy of this 
state with respect to the non interacting ground state is 
(AH) = (H ) + (V), with 
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where the sums on q and k are implicitly limited to q < 
Uf and k > kp. As we will check later (see below, eqn. 
©), <t>k,q ~ 1/fc 2 for large momenta, in order to satisfy 
the short range behavior 1 jr of the pair wave function 
in real space. This means that most of the sums on k 
diverge for k — > oo. This singular behavior is regularized 
by the renormalization of the coupling constant using the 
Lippman-Schwinger formula, thus yielding a vanishing 
qb- As a consequence, since the third sum in (V) is 
convergent, it gives a zero contribution to the final energy 
when multiplied by g#, and can therefore be omitted in 
the rest of the calculation. 

The minimization of (H) with respect to 0o and <j>k,q is 
straightforward and yields the following set of equations 
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where E is the Lagrange multiplier associated to the nor- 
malization of \ip), and can also be identified with the trial 
energy. These equations can be solved self consistently 
by introducing an auxiliary function x(l) — 0o+Efc 0fc-q 
and we obtain 
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After a straightforward calculation, this yields 
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where we got rid of the bare coupling constant gs by us- 
ing the Lippman-Schwinger equation ([2]). This equation 
can be solved numerically and yields E = —0.3h 2 k F /m, 
i.e. r\{i < —0.6 [26j. Note that the same analytical result 
was obtained independently in (27j . 
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FIG. 2: Comparison between experimental data of ref. [14I ] 
(squares) and eq. Q (full line). R a , Rp and Ri are re- 
spectively the radii of the superfluid, minority and major- 
ity components. The condensate vanishes for R a — at 
q = R}/Rl ~ 0.32. 



TRAPPED SYSTEM AND COMPARISON WITH 
EXPERIMENTS 

In the rest of the paper, we would like to show how 
experimental data from ref. [3] permits to improve the 
determination of the parameters T) a ,/3 ■ In this pursuit, 
we use the Local Density Approximation (LDA) to cal- 
culate the density profile of the cloud in a harmonic trap, 
which for simplicity we assume is isotropic. In this case, 
the chemical potentials of each species depends on 
position according to the law fiipfr) = n\ 2 — muj 2 r 2 /2, 
where uj is the trap frequency. 

Using this assumption, the transition between the var- 
ious phases will happen at radii R a and Rp given by 
^(Ra.js) / ^i(Ra,/3) = Va,/3- When these two equations 
are associated with the condition giving the radius R\ of 
the majority component, fii(Ri) = 0, we can eliminate 
both /LtJ and ^ from the equations, yielding the following 
close formula relating the radii R a , Rp and R\ 



Here, q 
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(4) 



of the R 2 p/R\ at which R a vanishes 



r)p) corresponds to the value 
i.e. at which the 

superfluid fraction disappears. In Fig. [21 we compare 
the prediction of eq. 2] with the experimental finding 
of ref. taking q = 0.32 to match the superfluidity 

thresholds. We see that close to the threshold, the agree- 
ment between the two graphs is quite good. However, 
they depart from each other for Rp/Ri > 0.7, corre- 
sponding to low population imbalance. One explanation 
for this discrepancy might involve finite temperature ef- 
fects. Indeed, it was already noted in Fig. 4 of ref. [l2j 
that, although the superfluid fraction was very sensitive 
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to temperature at small imbalances, the value of the crit- 
ical population imbalance was more robust. 

The superfluid phase disappears when R a vanishes. 
From cq. [H we see this happens for a ratio R 2 ^/R\ = q. 
As seen in Fig. [3J q can be extracted from the experimen- 
tal data of ref. [14|, which yield q ~ 0.32 and therefore 
constrain the possible values of rj a and rjp. Indeed, using 
this determination of q, as well as the rough upper and 
lower values for r] aj p, one obtains 
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0.10 < T] a < - 

-0.62 < r\p < 



0.088 
-0.60 



(5) 
(6) 



These bounds can be compared to the values deduced 
by BCS theory, predicting r\ a ~ 0.1 and r\p = 0. Our 
calculation excludes these values and explains why the 
width of the mixed normal state predicted by BCS theory 
is much narrower than observed in experiments. 



CONCLUSION 



In conclusion, we have presented an analysis of the ex- 



perimental data of ref. 14j providing stringent bounds 



on the values of the thresholds for quantum phase tran- 
sitions in uniform unbalanced fermi gases. Since they 
were obtained using minimal assumptions (mainly zero 
temperature and LDA), these bounds are fairly robust. 
In particular, they do not depend precisely on the su- 
perfluid nature of the intermediate phase. Our results 
suggest interesting follow-ups. First, the full understand- 
ing of the system, and in particular of the density profile 
of the cloud, requires the knowledge of the state equa- 
tion of the intermediate phase, whose exact nature then 
needs to be clarified. Second, the comparison with the 
data of ref. [l]| suggests an intriguing issue. Indeed, al- 
though the superfluidity threshold was not directly mea- 
sured in this paper, the parameter q can be inferred from 
the critical imbalance 0.7 measured by MIT. Rice's ex- 
perimental data yield at this value q ~ 0.16. Not only 
is this value very far from the one obtained here from 
the analysis of MIT's experiments, but it also contra- 
dicts the theoretical bounds rj a > —0.10 and np < —0.60 
which imply q > 0.31. As suggested in ref. [23j, this dis- 
crepancy may arise from surface tension effects provoked 
by the strong anisotropy of Rice's trap. Another inter- 
pretation might be the onset of the intermediate phase 
due to finite temperature effects, as suggested by some 
mean-field scenarios. Finally, the r\ a .p parameters can 
be evaluated experimentally using the value of q associ- 
ated with the measurement of the density discontinuity 
An^a at r = R a , given by Ani/An-2 = —rj a 24]. The 
preliminary data presented in fig. 2.b of ref. [2a ] suggest 
that the discontinuity Ani is very weak, hence indicat- 
ing a small value of rj a , in agreement with the bounds 
obtained here. 
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